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Summary 

This report describes the study of exhaust gas emissions from 
a staged Rich Burn/Quick-Quench Mix/Lean Burn (RQL) 
combustor under test conditions encountered in High Speed 
Research (HSR) engines. The combustor was modeled at con- 
ditions corresponding to different engine power settings, and 
the effect of primary dilution airflow split on emissions, flow 
field, flame size and shape, and combustion intensity, as well 
as mixing, was investigated. A mathematical model was devel- 
oped from a two-equation model of turbulence, a quasi-global 
kinetics mechanism for the oxidation of propane, and the 
Zeldovich mechanism for nitric oxide formation. A mixing- 
controlled combustion model was used to account for turbulent 
mixing effects on the chemical reaction rate. This model 
assumes that the chemical reaction rate is much faster than 
the turbulent mixing rate. The work described in this report is 
a part of the NASA High Speed Research Low NO x program. 

Introduction 

Future gas turbine engines are expected to operate at higher 
pressure ratios, inlet-air temperatures, and fuel-air ratios. At 
these conditions, conventional combustors would produce high 
levels of oxides of nitrogen (NO x ). Nitric oxides are divided 
into two classes according to their mechanism of formation: 
(1) “thermal NOx,” which is primarily nitric oxide (NO) 
with a small percentage of nitrogen dioxide (N0 2 ), is formed 
from the oxidation of diatomic nitrogen in the combustion air 
and fuel; and (2) “organic NO x ” is formed by oxidation of 
organically bound nitrogen in the fuel (fuel-bound nitrogen 
(FBN)). The main factors governing the thermal NO x 
formation rate are flame temperature and residence time; 
therefore, these factors must be lowered in order to reduce 
NO x formation. Organic NO x formation is less well under- 
stood than thermal NO x formation; however, we do know 


that measures which reduce flame temperature and thereby 
abate thermal NO x have little effect, positive or even adverse, 
an organic NO x . In this study, only thermal NO x is considered. 

Two promising concepts for reducing NO x emissions of 
advanced engines are the lean, premixed prevaporized (LPP) 
and the Rich Burn/Quick-Quench Mix/Lean Burn (RQL) 
combustors. The RQL concept (refs. 1 to 5) incorporates 
staged burning. Combustion is initiated in a fuel-rich zone with 
an equivalence ratio in the range of 1.2 to 1.8. The hydro- 
carbon reactions proceed rapidly, thereby depleting the available 
oxygen and inhibiting NO x formation. Additional air is then 
added, and combustion is completed under excess-air condi- 
tions. Combustion takes place in the rich and lean zones at 
temperatures well below the stoichiometric flame temperature, 
so NO x production is minimized. 

The lean, premixed prevaporized (LPP) combustor is another 
method of flame-temperature control, in addition to using 
variable geometry features to control the stoichiometry of the 
primary combustion zone (refs. 6 to 8). In this approach, 
increasing the amount of primary zone air decreases NO x 
emissions by reducing both flame temperature and residence 
time. 

In this study we consider a two-dimensional, confined 
swirling flow, staged RQL combustor model and study the 
combustor characteristics for a range of operating conditions. 
In the following sections the problem is formulated, the 
combustion models and numerical procedure are described, 
and the results are presented. 

Model Description 

The combustor consists of a primary and a secondary zone. 
The fuel (propane) and primary air are fed to the combustor 
via two confined, coaxial, swirling jets, with propane in the 
inner jet and air in the outer jet. Additional quick-quench air 
is fed to the secondary zone through the holes located on the 
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(a) Schematic of combustor model (dimensions are in meters). 

(b) Dilution jets velocity prolfile. 

Figure 1— Staged combustor and possible orientation of secondary swirling air jets. 


combustor liner. A schematic of the staged combustor, together 
with possible orientation of the secondary swirling air jets, 
is shown in figure 1. 

Governing Equations 

The governing equations for a turbulent axisymmetric, 
compressible, steady, and swirling flow case are considered. 
The general equations representing the conservation of mass 
and momentum for two-dimensional steady flows can be 
written in cylindrical coordinates (ref. 9) as continuity. 
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and the azimuthal (w) momentum component. 
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where the following modeling of fluctuating quantities has been 
used: 
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denote fluctuating quantities. The values of the turbulence 
(7) constants are given in table I. 
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Thermal and Reacting Flows 

The stagnation enthalpy equation can be written as 
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where the laminar and turbulence Prandtl numbers have been 
assumed equal to one, and the mean flow stagnation enthalpy 
H is defined as 


Turbulence Modeling 

The effect of turbulence is included by means of a k - 8 model. 
The differential equation for turbulence kinetic energy k is 
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In the foregoing equations, r and z are the radial and axial 
coordinates measured from the combustor centerline and 
combustor inlet, respectively; p is the mean flow density; u, 
v, and w are the mean flow axial, radial, and azimuthal velocity 
components, respectively; p is the pressure; P is the production 
rate of turbulence kinetic energy; and the effective viscosity 
/x e ff is the sum of the laminar (p) and turbulence (/z t ) viscosities. 
A scalar (isotropic) turbulence viscosity is assumed; C M , C|, 
C 2 , and o z are turbulence model constants, and the primes 


where Y. is the concentration of species j; hj is the enthalpy 
of formation of species j at the reference temperature 
T 0 = 298. 15 K; C p is the specific heat at constant pressure; 
T is the temperature; and the five chemical species considered 
are C 3 H 8 , 0 2 , N 2 , C0 2 , and H 2 0. 

The combustion of propane is modeled by using the one- 
step global mechanism: 


C 3 H 8 + 5(0 2 + 3.76N 2 ) - 3C0 2 + 4H 2 + 18.8N 2 (18) 

The oxidation of fuel is controlled by the chemical reaction 
rate. This reaction rate is controlled either by the mixing of 
eddies (mixing-controlled) that contain reactants and products 
or by an Arrhenius rate equation (chemical kinetics controlled). 
In this investigation, the mixing-controlled model is used. In 
the mixing-controlled model, the influence of turbulence on 
the reaction rate is taken into account by assuming that the 
generation and the dissipation of the time-average of the square 
of the fuel concentration fluctuations are equal: 


Cg t pk u2 t 



where T, is the mean fuel (propane) mass fraction, y,- is the 
fluctuating mass fraction of fuel, and Cg\ and Cg 2 are 
constants. 

The mean reaction rate of the fuel /?, can be written as 
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TABLE I. -VALUES OF MODEL CONSTANTS 


Turbulence model 

Mixing-controlled model 

c. 

C, 

c 2 

Cg\ 


Q?3 

C 

0.09 

1.44 

1.0 

2.70 

0.134 

0.00456 

0.02 


where 



The values of the mixing-controlled reaction rate model 
constants are shown in table I. 


where R is the universal gas constant and W denotes molecular 
weight. 

Boundary Conditions 

At the combustor inlet the How was assumed to have a zero 
mean radial velocity, and uniform axial velocity and pressure. 
The azimuthal velocity profile of the inner (fuel) jet was 
simulated by a solid body rotation, whereas that of the outer 
(inlet-air) jet was modeled by a free vortex. The temperatures 
of the fuel and inlet-air jets were assumed to be uniform and 
equal to 300 and 1000 °F, respectively. The inlet profiles of 
turbulence kinetic energy taken from references 10 and 1 1 are 
as follows: 


Species Equation 

If the laminar and turbulence Schmidt numbers for all 
species are assumed equal to one, the following equations are 
obtained for the time-averaged mass fractions: 
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where 
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5 is the reaction stoichiometric coefficient; Y 0 , T n , and T p are 
the mass fractions of the oxidizer, nitrogen, and combustion 
products, respectively; and Y sz is the mass fraction of the 
Shvab-Zeldovich variable used in the Zeldovich mechanism 
for nitric oxide formation. 


Equation of State 

The density can be calculated from the equation of state: 
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where R, and R () are the inner radii of the fuel and inlet-air 
pipes, respectively; R s is the outer radius of the fuel pipe; and 
Uj and u () are the mean axial velocities of the fuel and inlet- 
air jets, respectively. The combustor geometry and dimensions 
are shown in figure 1(a). 

At the combustor centerline (symmetry axis) both the radial 
and azimuthal velocity components are zero, whereas the radial 
derivatives of the other How quantities are zero. At the 
combustor outlet, the flow is assumed to leave with zero axial 
derivatives. At solid walls the normal velocity component to 
the walls was set to zero, and the law-of-the-wall was used 
to calculate the tangential velocity near the walls. The solid- 
wall temperature was assumed equal to 350 K. The Reynolds 
analogy was used to calculate the heat transfer loss through 
the solid walls. 

The boundary conditions of the quick-quench mixing jet and 
the dilution jets can be summarized as follows: 

u } = V sin a: cos 0 

Vj = V sin a sin 


vv’j = V cos a 

where u r vj, and are the velocity components of the jets 
and V, a, and (3 are shown in figure 1(b). Both of the jets’ 
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mass How rates and pressures are prescribed. The inlet profiles 
of the jets’ turbulence kinetic energy are similar to that of the 
combustor’s inlet fuel jet. 

Numerical Method 

The numerical procedure employed in the calculations herein 
is of the control volume variety (ref. 12). It uses a staggered 
grid and an implicit finite-difference algorithm. Calculations 
were made by a hybrid numerical technique that uses control 
differences for the convection terms whenever the absolute 
value of the mesh Reynolds number is smaller than two, and 
uses upwind differences when the absolute value is greater than 
two. The calculations were performed in the following manner. 
First, a value for an initial pressure field was assumed and 
employed to calculate the corresponding velocity field. Then 
a Poisson equation for the pressure field was formulated and 
taken to be equal to the sum of the assumed field value plus 
a local pressure perturbation. The solution of the Poisson 
equation yields the new pressure, which is generally different 
from the guessed one. Therefore, the momentum and Poisson 
equations were solved iteratively until a prescribed 
convergence criterion (error tolerance) was satisfied. A grid 
consisting of 32 by 32 points was used in the calculations. 

Results and Discussion 

Calculated results are presented here for the Rich 
Burn/Quick-Quench Mix/Lean burn combustor. The 13 
different How conditions that were studied are listed in table II. 
In figures 2 to 8 the profiles of the velocity vectors, stream 
function density, turbulence intensity (/ = \f2k/3/u i ), gas 
temperature, NO x distribution, and fuel, oxygen, nitrogen, 
carbon dioxide, and water distributions for <F RZ (rich-zone 
equivalence ratio) = 1.65 (case 8) are presented. Figure 2 


presents the mean velocity vectors showing magnitude and 
direction for a counter-swirl condition, that is, the fuel and 
inlet-air jets are counter-swirling. The inlet How is lifted by 
centrifugal forces and creates a recirculation zone at the 
combustor centerline. The recirculation zone is characterized 
by a steep axial and azimuthal velocity gradient. A negative 
velocity, which appears at the location where the recirculation 
zone starts rolling up, disappears as the recirculation zone 
approaches the combustor inlet. The inlet-air flow shows a 
more uniform axial velocity profile along the recirculation zone 
and has a very steep velocity gradient at the shear layer regions 
of the incoming jets. The radial distribution of the velocity 
flattens out and becomes uniform as the velocity differs radially 
and axially downstream of the recirculation zone along the 
combustor. Figure 2 also shows that the quick-quench mixing 


TABLE II.-RQL COMBUSTOR OPERATING 
CONDITIONS 

(Inlet temperature, 1000 K; inlet pressure, 2.0 atm.] 


Case 

number 

Rich-zone 

equivalence 

ratio, 

0RZ 

Quick-quench 

zone 

equivalence 

ratio 

Lean-zone 

equivalence 

ratio 

1 

1.09 

0.81 

0.79 

2 

1.19 

.83 

.81 

3 

1.21 

.78 

.76 

4 

1.28 

.84 

.82 

5 

1.40 

.95 

.93 

6 

1.43 

.78 

.76 

7 

1.50 

.72 

.70 

8 

1.65 

.79 

.77 

9 

1.80 

.86 

.84 

10 

2.09 

.89 

.87 

11 

1.50 

.77 

.75 

12 

1.50 

.80 

.78 

13 

1.50 

.82 

.80 
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0 
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Figure 2.— Mean velocity vector profile for a rich-zone equivalence ratio of 1.65. 
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jet penetrates radially into two-thirds of the combustor radius, 
and that a portion of the quick-quench mixing jet recirculates 
upstream and creates a pair of lopsided vorticies sandwiching 
the jet. The dilution jet penetration is much less pronounced 
since the dilution jet-to-mainstream momentum flux ratio is 
lower than the quick-quench jet-to-mainstream momentum flux 
ratio. Comparison with the stream function density profile 
shown in figure 3 confirms the flow structure and dynamics 
observed from the velocity profile. 

Figure 4 presents the turbulence intensity profile along the 
combustor. The turbulence intensity at the shear layer regions 
is much smaller than that in the recirculation zone. The high 
turbulence intensity observed in the recirculation zone and the 
lean-zone region downstream of the quick-quench mixing jet 


can be attributed to the high turbulence viscosity in these 
regions. The large radial gradient of the azimuthal velocity 
component in the recirculation zone results in a high production 
and dissipation rate of turbulence kinetic energy in this region. 
The high diffusion of axial and azimuthal velocity components 
in the lean-zone region downstream of the quick-quench 
mixing jet accounts for the large generation and dissipation 
rates of turbulence kinetic energy in this region. Turbulence 
is generated in regions of steep velocity gradients, for example, 
at the shear layer regions, since the turbulence kinetic energy 
produced is proportional to the square of the velocity gradients. 

Figure 5 shows the mean temperature profile of the gas along 
the combustor. Note that the temperature is high in regions 
with high turbulence intensity. This is due to better mixing 



Stream function density 



Figure 3.— Stream function density profile for a rich-zone equivalence ratio of 1.65. 



1.27 

Turbulence intensity 



Figure 4.— Turbulence intensity profile for a rich-zone equivalence ratio of 1.65. 
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and the higher reaction rate that occurs in these regions. Most 
of the reactions take place at the rear of the high-pressure core 
of the recirculation zone and in the fuel-lean region 
downstream of the quick-quench mixing jet. Quenching takes 
place in the quick-quench mixer region because of the high- 
flow rate, quick-quench mixing air injected through the single 
hole on the two-dimensional combustor wall. An overall 
equivalence ratio of 0.5 is thus achieved. This results in very 
low combustion efficiency in this region. 

The NO x distribution along the combustor is shown in 
figure 6; it is highest in the lean-zone region, which would 
be expected since NO x formation is proportional to flame 
temperature, residence time, and the availability of oxygen. 
The formation rate of NO x is much smaller than the oxidation 
rates of hydrocarbon fuel and CO. 

Figure 7 shows fuel vapor, oxygen, and nitrogen distri- 
bution. The fuel-rich region located at the front end of the 


recirculation zone accounts for a temperature that is lower in 
this region than at the rear end of the recirculation zone. 
Comparison of the fuel, oxygen, and nitrogen distributions 
confirms the observations made from the temperature and 
NO x profiles. Figure 8 presents the carbon dioxide and water 
distributions along the combustor. The distribution of these 
combustion products agrees with the predicted flame 
temperature and fuel and oxygen profiles. 

In figure 9 we summarize the data on NO x concentration 
as a function of the rich-zone equivalence ratio. These data 
indicate that as the rich-zone equivalence ratio increases from 
1.19 to 2.0, the NO x concentration decreases. 

Figure 10 shows the radial profile of combustor exit 
temperatures for cases 11 to 13. The maximum temperature 
is highest for case 13. The combustor liner heat transfer 
profiles are shown in figure 1 1 for cases 11 to 13. 
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Figure 5.— Temperature distribution for a rich-zone equivalence ratio of 1.40. 



NO x mass fraction 



Figure 6. — NO x mass fraction distribution for a rich-zone equivalence ratio of 1.40. 
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(a) Fuel vapor distribution. 

(b) Oxygen distribution. 

(c) Nitrogen distribution. 

Figure 7.— Fuel vapor, oxygen, and nitrogen distributions for a rich-zone equivalence ratio of 1.40. 
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Carbon dioxide 




Low High 

Water 



(a) Carbon dioxide distribution. 

(b) Water distribution. 

Figure 8.— Carbon dioxide and water distributions along the combustor for a rich-zone equivalence ratio of 1.40. 
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Rich-zone equivalence ratio, (j)^ 

Figure 9.— NO x concentration as a function of rich-zone equivalence ratio 
for RQL combustor. 



Outlet temperature, K 

Figure 10.— Typical outlet temperature profiles for RQL combustor. 


Concluding Remarks 

A mathematical model and a computer program were 
developed to study the How field, combustion, and exhaust 
gas emissions of a low NO x , staged combustor representative 
of the Rich Burn/Quick-Quench Mix/Lean Burn concept. The 
two-dimensional axisymmetric combustor consisted of a 
primary and a secondary zone. 

The model solved the two-dimensional equations of 
continuity, momentum, energy, and chemical species as well 



Axial distance, m 


Figure 1 1.— Typical local liner heat transfer profiles for RQL combustors. 


as a two-equation model of turbulence. Discrete equations were 
solved by the semi-implicit method for pressure-linked 
equations (SIMPLE). Turbulent boundary layers were resolved 
with law-of-the-wall equations, and wall heat fluxes were 
calculated from the Reynolds analogy between momentum and 
heat transfer. The oxidation of the fuel was modeled by using 
a quasi-global kinetic mechanism, and this reaction rate was 
modeled by a mixing-controlled model. The amount of nitric 
oxide formation was computed for a variety of operating 
conditions by using the Zeldovich mechanism. 

The calculations indicate that high turbulence and mixing 
yield high reaction rates and high flame temperatures in 
addition to good combustion efficiency. The predicted NO x 
formation is a function of flame temperature, residence time, 
and oxygen availability. Our findings, however, show that 
quick-quench mixing and dilution jets in gas turbine com- 
bustors constitute a fully three-dimensional problem that must 
be simulated by using a three-dimensional cylindrical 
coordinate system. The two-dimensional simulation of the 
quick-quench mixing and dilution jets in this study results in 
local quenching of the reaction in the quick-quench region, 
which is due to the high flow rate of air injected through the 
two-dimensional hole on the combustor liner. Future work will 
focus on developing a three-dimensional version of the present 
technique and extending the mixing-controlled model to 
include more chemical species and chemical reactions. 
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